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The multifaceted role of the density dependent nuclear symmetry energy in the nuclear astro- 
physics involving neutron stars is highlighted. Efforts toward a model independent determination 
of the dense matter equation state through a deconstruction of the neutron star structure equa- 
tion utilizing the masses and radii of several individual neutron stars are described. The need for 
observational data of both measurements for the same star is stressed. 
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1. Scope of This Talk 

In this talk, I will 

1. summarize some theoretical advances made in establishing correlations between observables 
of nuclei and neutron stars, 

2. indicate how laboratory experiments, e.g., at JLab and Rare Isotope Accelerators (RIA's), 
can help to unravel the composition and structure of a neutron star, 

3. point out how the observed masses (M's) and radii (7?'s) of several individual neutron stars 
can be used to construct the model independent equation of state (EOS) of dense matter 
(through an inversion process that I term as "Deconstructing a Neutron Star"), and 

4. indicate the key neutron star observations that are necessary to take leaps in our understand- 
ing of the nature of strong interactions at supra-nuclear densities. 



2. The Pervasive Role of the Nuclear Symmetry Energy 



The nuclear (a)symmetry energy 

I d^e{n,x) nn-Up 
E,ym[n,5) = - — ' ^ = — T — = l-2x, (2.1) 

represents the energy cost required to create an isospin asymmetry 8 in nucleonic matter. Above, 
Hn and Hp are the neutron and proton densities, x is the proton fraction in matter, and e{n,5) is the 
energy density of isospin asymmetric matter. Figure 1 from Ref. [|l|] shows the energy per particle 
of bulk nuclear (x = 1 /2) and neutron matter (x = 0) versus the total baryon density n. 




Figure 1: Results shown are for the microscopic potential model calculations of Akmal and Pandharipande 
(APR) |[] and simpler nonrelativistic and relativistic effective model fits (NRAPR and RAPR) [[|]. 
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The physical properties of nuclei, such as their masses, neutron and proton density distribu- 
tions (including their mean radii), collective excitations, fission properties, matter and momentum 
flows in high energy heavy-ion collisions, etc. all depend on the isospin structure of the strong 
interactions between nucleons (i.e., nn and pp interactions versus np interactions). The energy 
jl = jjLn — jjLp = AEsymd, whcrc jXn and jXp are the neutron and proton chemical potentials, respec- 
tively, is crucial in determining reaction rates involving electrons and neutrinos, particle abun- 
dances, etc., in astrophysical contexts such as supernova dynamics, proto-neutron star evolution, 
the r— process, the long-term cooling of neutron stars, and the structure of cold-catalyzed neutron 
stars (i.e. their masses, radii and crustal extent), etc. The pervasive role of the isospin dependence 
of strong interactions in nuclear processes in the laboratory and the cosmos is sketched in Fig. ^ 
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Figure 2: The pervasive role of the nuclear symmetry energy. Figure taken from Ref. m 



3. Connections to Neutron Star Structure 

The structure of neutron stars is determined by the energy and pressure of charge-neutral beta- 
stable matter, the nucleonic components of which can be written as 

E{n,x) ^ E{n,0.5)+E,y,n{n) {l-2xf + ---, 

Pin,x) ^n^[E'in,0.5)+E',,^{l-2x)^]+---, (3.1) 
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where the primes above refer to derivatives with respect to density. Leptonic (electrons and muons) 



contributions must be added to Eq. (|3JJ) for the total energy and pressure. As demonstrated in 
Ref. [Q], the radius (but, not necessarily the mass) of a nucleonic neutron star is chiefly determined 
by the pressure in the vicinity of the nuclear equilibrium density, uq ~ 0. 16 ± 0.01 fm^^, for which 



the contribution of the first term to the pressure in Eq. ( [3.1| ) is very small. This feature of a neu- 
tron star's structure highlights the importance of the density dependence of the nuclear symmetry 
energy near hq. Figure |3| shows the mass-radius relationships for select EOS models including the 
cases in which extreme softening due to the presence of exotica such as hyperons, quarks or Bose 
condensates occurs in the cores of stars. 
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Figure 3: Mass versus radius curves from Ref. [|]|. EOS symbols are as in Ref. 



In Fig. 0, the observed neutron star masses are shown. The simple and error weighted means 
for neutron stars in x-ray binaries are 1.55 Mq and 1.37 Mq, respectively. The corresponding 
numbers for double neutron star binaries are 1.32 M0 and 1.41 Mq, whereas for white dwarf - 
neutron star binaries, the numbers are 1.60 Mq and 1.33 Mq. It is unfortunate that in those cases 
(e.g., binary radio pulsars) in which the masses are very accurately known, the corresponding 
numbers for radii are unknown. It is hoped that with improvements in timing analysis, the moment 
of inertia of the neutron star in the double-pulsar system J0737-3039 will become available [^. As 
will be shown later, the masses and radii of several individual neutron stars can uniquely pin down 
the EOS of dense matter. For prospects of obtaining this much needed information, see the article 
by Bob Rutledge in these proceedings. For prognosis of EOS determination, see Ref. [^. 
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Figure 4: Observed neutron star masses. Figure courtesy of Lattimer. 



4. The Role of Isospin Interactions in Nuclei and Neutron Stars 

A few of the recently discovered empirical relationships that underscore the role of isospin 
interactions in nuclei and neutron stars are highlighted here. 

The neutron skin thickness in nuclei and the sub-nuclear pressure of pure neutron matter 

Typel and Brown M, have noted that model calculations of the difference between neutron 
and proton rms radii 5R = {r^) — {rf,) are linearly correlated with the pressure of pure neu- 
tron matter at a density below no characteristic of the mean density in the nuclear surface (e.g., 
0.1 fm^-'). The density dependence of the symmetry energy controls the so-called neutron skin 
thickness 5R in a heavy, neutron-rich nucleus. Explicitly, 5R is proportional to a specific average 
of [1 — Esym{no) / Esym{n)] in the nuclear surface, see Refs. | [T0| , |I|]. 

The left panel of Fig. || shows the Typel-Brown correlation extended to neutron-star matter 
in Ref. [|l|] using both potential and field-theoretical models. The proposal [ pTj , p^ ] at JLab to 
measure the neutron rms radius of ^o^p^, through parity-violating electron scattering experiments 
to 1% accuracy (the rms charge radii of nuclei are known to better accuracy) can help to provide a 
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Figure 5: Left panel: The neutron skin thickness 5R of finite nuclei versus the pressure of /3 -equilibrated 
neutron star matter at a density of 0.1 fm^^^. Right panel: Calculated neutron skin thicknesses 5R of nuclei 
versus the radii of 1.4M(;) stars. Trends for stars up to the maximum mass are similar. The solid lines 
represents a least square fit. Trends for stars up to the maximum mass are similar [[l|]. Figure taken from [|l|]. 



calibration pressure of neutron-star matter at sub-nuclear densities. The knowledge of neutron rms 
radii for neutron-rich nuclei of varying masses will be of great help in establishing the uncertain 
density dependence of isospin interactions at sub-nuclear densities. The right panel of Fig. 5 shows 
that the correlation between 5R and R14 is not very sharp. Accurate predictions of neutron star 
radii require knowledge of the equation of state at supra-nuclear densities. 

The neutron skin thickness in nuclei and the neutron star radius 

Horowitz and Piekarewicz have pointed out that models that yield smaller neutron skins in 
heavy nuclei tend to yield smaller neutron star radii. These authors, along with others [14, I]], have 



stressed the need for accurate measurements of the neutron skin thicknesses of nuclei. Although 
the neutron star radii are determined at supra-nuclear densities, the calibration provided at sub- 
nuclear densities cannot be ignored. Even with varying stiffness at high density, several works 
have confirmed the trend shown in the right panel of Fig. ^ The Horowitz-Piekarewicz correlation 
is thus very suggestive inasmuch as radii of stars up to the maximum mass exhibit similar trends. It 
would be fruitful to test this correlation from accurate radius measurements of neutron stars whose 
corresponding masses are also well determined. 

The neutron star radius R and the pressure P of neutron-star matter 

Lattimer and Prakash [Q] found that the quantity RP^^I^ is approximately constant, for a given 
neutron star mass, for a wide variety of EOS's when the pressure P of beta-equilibrated neutron- 
star matter is evaluated at a density in the range mq to 2«o, where hq denotes equilibrium nuclear 
matter density. Since the pressure of nearly pure neutron matter (a good approximation to neutron 
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star matter) near hq is approximately given by n E[y,^ (see Eq. (3.1)), the density dependence of 
the symmetty energy just above no will be a critical factor in determining the neutron star radius. 
Figure^ shows results for \AMq stars. Similar trends persist for stars up to the maximum mass star 
as shown in Ref. [Q]. This correlation, coupled with accurate measurements of neutron star radii, 
can delimit the very uncertain (up to a factor of 5!) pressure of neutron-star matter (see Ref. [Q]) at 
supra-nuclear densities. 
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Figure 6: The quantity RP^" as a function of the stellar radius R, for pressures P determined at 3/2, 2 
and 3 times equilibrium nuclear matter density. For each density, the best-fit value for the exponent a is as 
indicated. Circled stars indicate the results obtained with the APR equation of state. These results are for 
1.4M0 stars; results for stars up to the maximum mass star are similar. Figure taken from Ref. 



5. Deconstructing a Neutron Star 



This section contains a report of on-going work of the speaker with Sergey Postnikov and 



James Lattimer [15]. The question being posed here is "Can observations of the masses and radii 



of several (say 5 to 7) individual stars uniquely determine the dense matter EOS?" In short, the 



answer is a resounding YES ! This assertion follows from the work of Lindblom [ |16| ] who exploited 
the one-to-one correspondence between an EOS and the M — R curve generated through the use of 
the Tolman-Oppenheimer-Volkov (TOV) [ 17 , equations of stellar structure. Given the diverse 
theoretical predictions, it is of great interest to enquire whether or not observations of masses and 
radii of several neutron stars can be used to uniquely pin down the poorly known high density 
EOS. The extent to which inherent errors in measurements affect the extraction of the EOS (not 
addressed in Lindblom's work) will become an important issue. 



Recently launched campaigns [19] to accurately measure masses and radii of a large number 
of neutron stars give impetus to investigate this inversion procedure in more depth than originally 
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undertaken by Lindblom. The extension of such an inversion procedure to include other observ- 
ables such as the surface redshifts, binding energies, moments of inertia, etc. has been undertaken 
by Postnikov, Prakash and Lattimer incorporating the inherent errors involved in observations. 

We begin by recasting the TOV equations of stellar structure using the variable h defined 
through dh = dp/{p + p (p)), where p is the pressure and p {p) is the mass-energy density. Adopt- 
ing the units G = 1 and c = 1, the result is 

dr^ 9 r — 2m dm ^ r — 2m 

-fr = -^r^ — ^) — = -4nPp — ^, (5.1) 

dh m + 4Tlr^p dn m + ^nr^p 

where p{h) and p{h), which serve as input, contain the EOS. The advantages of this reformulation 
(somewhat different than that used by Lindblom) are that the enclosed mass m and radius r are now 
dependent (on h and thus the EOS) variables, and h is finite both at the center and surface of the 
star. Furthermore, the quantities r^{h) and m{h) admit tractable Taylon expansions about he that 
feature the finite quatities pc, Pc and {dp/dh)c, where the subscript c denotes the star's center. 

The procedure is to begin with a known EOS up to a certain density, take small increments in 
mass and radius, and adopt an iterative scheme based on Eq. ( ^ ) to reach the new known mass 
and radius. In the absence of observational numbers for individual stars, a test of the above scheme 
can be performed by assuming that the masses and radii are those that result from a model EOS. 

We have tested our scheme using an EOS consisting of two polytropes whose indices are 
appropriately chosen to mimic a realistic EOS and obtained results that reproduce the exact EOS 
to hundredths of percent accuracy. The left panel of Fig. ^ shows our results for this test case. 
Methods to yield even better accuracy using this modified Lindblom procedure are being devised. 

An alternative procedure is to solve the reformulated TOV equations from the center to the 
surface with the assumed form of the EOS using a Newton-Raphson scheme to obtain the known 
mass and radius. The results from this approach is shown in the right panel of Fig. ^. This method 
also yields results to hundredths of percent accuracy. We now have two methods that yield similarly 
accurate results. 




Figure 7: Deconstructing a neutron star with a two-power polytropic EOS. 



We are now in a position to test our two schemes using physically motivated model EOS's. As 
examples of realistic EOS's, we will begin with the EOS's of Ref. pO| ] . These EOS's allow for 
variations in the uncertain high-density behavior of the nuclear symmetry energy and mimic trends 
found in microscopic calculations of high-density matter. The EOS labeled PALS 1 is characterized 
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by the nuclear matter compression modulus of ^ = 240 MeV and a symmetry energy function for 
potential interactions: F{u) = u, where u = «/«o and uq = 0.16 fm^^. Including contributions from 
kinetic energy, the symmetry energy at «o is 30 MeV. The results in Fig. || show the extent to which 
the EOS is reconstructed for the case in which the masses and radii are assumed to be those that 
result from the EOS of PAL31 [|^]. The two iterative schemes described above yield results of 
similar satisfactory accuracy. 




Figure 8: Deconstructing a neutron star with a physically motivated nucleonic EOS. 



Work is in progress to determine the extent to which "Deconstruction" works accurately for a 
variety of EOS's including those that permit extreme softening at high density. Several exercises 
can be envisaged: (i) To what extent can one determine if only one mass and radius are very well 
determined? (ii) What if several masses and radii are known, albeit with unsatisfactory error bars? 
(iii) To what extent can we determine the EOS given the inherent error bars in measurements?. 
Answers to these questions and many others will be published soon [15|. In addition to its inherent 
worth, this investigation will aid astronomers to better plan and secure observations in the future. 



6. Conclusions 

It is clear that 



1. the masses and radii of several (say 5 to 7) of individual neutron stars can pin down (through 
deconstruction) the equation of state of neutron-star matter and shed light on the density 
dependence of the symmetry energy, i.e., the poorly known isospin dependence of strong 
interactions, and 

2. precise laboratory experiments, particularly those involving neutron-rich nuclei, are sorely 
needed to pin down the near-nuclear aspects of the symmetry energy which determines the 
masses, neutron skin thicknesses, collective excitations, etc., of nuclei. 



It is also clear that continuing mutual interactions of astronomers, laboratory experimenters 
and theorists in the fields of astrophysics and strong interaction physics are necessary to fulfill the 
objective of understanding the physics of compact objects in which the ultimate energy density of 
observable cold-baryonic matter is realized [21]. 
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